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Abstract 

A large class of multidimensional nonlinear Schrodinger equations admit localized nonradial 
standing wave solutions that carry nonzero intrinsic angular momentum. Here we provide 
evidence that certain of these spinning excitations are spectrally stable. We find such waves for 
equations in two space dimensions with focusing-defocusing nonlinearities, such as cubic-quintic. 
Spectrally stable waves resemble a vortex (non-localized solution with asymptotically constant 
amplitude) cut off at large radius by a kink layer that exponentially localizes the solution. 

For the evolution equations linearized about a localized spinning wave, we prove that un- 
stable eigenvalues are zeros of Evans functions for a finite set of ordinary differential equations. 
Numerical computations indicate that there exist spectrally stable standing waves having central 
vortex of any degree. 
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1 Introduction 

A long-term goal of the study of nonlinear waves is to sort out what kinds of structures and interac- 
tions are significant and robust. From this point of view, nonlinear Schrodinger equations are worthy 
of study because of their status as a canonical model for weakly nonlinear phenomena in many fields 
(including nonlinear optics as a prominent example), and as a simple model for more complicated 
gauge theories of mathematical physics. In more than one space dimension such equations admit 
a variety of solutions with interesting structure, including solitary waves (e.g., "light bullets") and 
vortices (kin to the magnetic vortices of the theory of superconductivity) . 
Here we consider nonlinear Schrodinger equations of the form 

— i dtu — Au = g(u) (1-1) 

where u:M. 2+1 — > C, and where g(u) = /i(|u| 2 )u for some real-valued C 1 function h. Our most 
significant results concern nonlinearities for which focusing (h! > 0) at small amplitude is overcome 
by defocusing (h! < 0) at larger amplitude. Prototypical is the cubic-quintic nonlinearity 

g(u) = \u\ 2 u- \u\ 4 u. (1.2) 
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Much of the theory developed in this paper also applies to equations of focusing type. 

We are interested in the stability of solitary- wave solutions with the symmetry of a vortex, having 
the special form 

u (x,t) =e wt+lm6 w{r) (1.3) 

in terms of polar coordinates (r, 9) for K 2 . Here w is a real constant (the standing wave frequency), 
m is an integer (the vortex degree), and w: [0, oo) — * R. We also refer to m as the spin of the 
wave, because the conserved quantity L that corresponds via Noether's principle to the rotational 
invariance of ( |l . l| ) is the (third component of the) angular momentum 

(—ix x Vu)u. 



For standing waves having the form (L3), L = 2-kvh J q 



The standing-wave amplitude profile w in ( |l.3[) must satisfy the differential equation 

I Ul 2 

w" + -w' T w + f(w) = 0, (1.4) 

where f{w) — g(w) — low. It is known that, under suitable conditions on /, for each integer m this 
equation has C 2 solutions w such that w(r) ^ exponentially as r — * oo. Iaia and Warchall fl5| 
proved that for fixed m, there is at least one such solution w for each prescribed number n of positive 



zeros. Each such solution gives rise to a localized standing wave solution of (LI). In turn, this gives 
rise to solitary-wave solutions traveling at any speed, via the Galilean boost transformation for the 
Schrodinger equation. 

Without loss of generality we assume h(0) — 0. Then the restriction of g to the real axis is a C 1 
function with <?'(0) = and g'(s) = 0(s 2 ) as s — > 0. For an exponentially localized solution to exist, 
it is necessary that /'(0) < 0, so the standing- wave frequency to is positive. It is also necessary that 
F(s) = J Q f(a) da be positive for some positive value of s attained by profile amplitudes w(r). 

In this paper we investigate the linear stability of such exponentially localized solitary waves. 



Linearizing (1.1) about such a solution yields a real system of the form 9 t $ = A<& where A is a 2 x 2 
matrix of second-order partial differential operators that has purely imaginary essential spectrum 
c C ss = ■ \t\ > lo}. Any eigenvalues of A not in the essential spectrum are discrete (isolated and 
of finite multiplicity) and occur with Hamiltonian symmetry: If A is an eigenvalue, then so are —A 
and A. 

We prove that the search for unstable eigenvalues of A reduces to the study of a finite number of 
eigenvalue problems for 2x2 systems of second-order ordinary differential equations. Eigenvalues 
for each of these systems correspond to zeros of an Evans function or Wronskian that we prove is a 
globally analytic function in the cut plane C CT = C \ a css . The real part of any unstable eigenvalue 
is bounded by a quantity determined from the profile shape. 

We use the argument principle and numerical computation to locate zeros of the Evans func- 
tions. Very often the operator A has unstable eigenvalues (meaning eigenvalues with positive real 
part), indicating that the associated solitary wave is linearly exponentially unstable with respect to 
perturbations of initial data. 

However, we have discovered spectrally stable waves (waves with no unstable eigenvalues) in a 
certain parameter regime for several nonlinearities of focusing-defocusing type, including the cubic- 



quintic (L2). All of these waves have radial profiles w(r) with no positive zeros. We have found such 
waves for all values of spin index satisfying \m\ < 5, with indications that they exist for arbitrary 
m . 

The spectrally stable spinning waves (m ^ 0) are found when to is close to but less than a critical 
value for which the two bounded regions between the graphs of y = g(x) and y = lu*x for x > 
have equal area. This condition on to* is equivalent to the condition that the potential function F 
with to = lo* has a zero- height local maximum or "hilltop" at a point a* > 0. For the cubic-quintic 
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case ojs,. = jjr = 0.1875 and spectrally stable waves are found for co in an interval uj a {rri) < u> < 
which shrinks as \m\ increases, with width roughly proportional to |to|~ 2 . 

A transition to instability occurs as the standing wave frequency u> decreases below lu ci . Un- 
stable eigenvalues appear, for a particular azimuthal mode, due to collisions of pairs of imaginary 
eigenvalues as in a Hamiltonian Hopf bifurcation. 

The waves develop a distinctive structure as lo approaches from below. (See Figs. |^ and 
U in section 2.) The wave amplitude w(r) increases from zero to remain approximately constant 
over a large r-interval, taking values near a*, near the hilltop of the nonlinear potential F, before 
rapidly returning exponentially to zero. Thus the waves have a structure consisting of a core region 
resembling a vortex (a non- localized solution with asymptotically constant amplitude, cf. [p2|). 
encapsulated or cut off (at an apparently arbitrarily large radius) by a circular kink region that 
serves to exponentially localize the solution. 

What are the salient features of the nonlinearity that might explain the existence of stable 
encapsulated vortices with this kind of structure? We believe that a basic understanding begins 
with the stability of spatially constant solutions 

u(x,t) = ae lut (1.5) 

where a > is constant and w — g(a)/a. This solution is linearly unstable if h'(a 2 ) > (the 
focusing case) and is linearly stable if h'(a 2 ) < (defocusing) ; see j2^j and Remark |3.2| below. 
Equivalently, stability is determined by the sign of g' '(a) — g(a) /a, a quantity that is easily visualized 
on the graph of g over the positive reals. The spatially constant solution ae M * is linearly unstable 
if g'(a) — g(a)/a > and linearly stable if g'(a) — g(a)/a < 0. 



A focusing-defocusing nonlinearity like the cubic-quintic (1.2) has two features that appear to 
be key (see Fig. |l|): 



(i) There is a number ao > such that the constant-state solution (1.5) is linearly unstable for 
< a < ao (h'(a 2 ) > 0) and is linearly stable for a > ao (h'(a 2 ) < 0). 

(ii) There is a number > for which the graphs of y — g(x) and y — uj*x enclose two regions of 
equal area, and so for which the potential F has zero- height hilltops at s = and at s = a* > 0. 

As we discuss in section 6, for any nonlinearity with these properties there is a linearly stable one- 



dimensional kink solution of (1.1) that connects the state a*e iw ** to the zero state. The value 



is the standing-wave frequency of the kink. In the limit as uj approaches and |m| — > ao, we will 
see that the structure of the encapsulated vortices consists of parts that approximate linearly stable 
solutions (constant states or kinks). In this way it becomes plausible that these waves are linearly 
stable. 

The nonlinear Schrodinger equation is of fundamental importance in nonlinear optics, where a 
cubic (Kerr) nonlinearity is most often assumed. A substantial body of work exists on so-called 
saturable media, for which the nonlinearity takes the form 

0(u) = Jh£L. (1.6) 
1 + \ur 



The cubic-quintic nonlinearity (1.2) can be obtained by truncating the Taylor expansion for this 
nonlinearity, but in fact there is a substantial difference between the saturable and cubic-quintic 
nonlinearities at large amplitude. This difference is reflected in the stability of spinning solitary 
waves. In the saturable case the spinning waves with |m| = 1 have been found to be unstable 



[|12|. Some spinning waves with \m\ = 1 for the cubic-quintic nonlinearity (1.2), however, have 
previously been found to be stable in direct numerical simulations by Quiroga-Teixeiro and Michinel 
p7f . We conjecture that the crucial difference is that the cubic-quintic is focusing-defocusing, but 
the saturable nonlinearity is only focusing — it does not admit any linearly stable nonzero constant 
states. 
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There is experimental evidence that the cubic-quintic model is in fact appropriate for certain 
optical materials (lj, [I?], [jl], [n|]. Instabilities of three-dimensional spinning waves for the cubic- 
quintic model have been investigated numerically in J9|, p0| . After this paper was first submitted, 
we learned of work of Towers et al. on the two-dimensional cubic-quintic model |u| . These authors 
numerically obtain stability results essentially similar to ours, with one discrepancy concerning the 
existence of a very weak instability for a range of values of u> in cases for |m| = 1,2,3 where we find 
stability. For related work on solitary waves with vortex symmetry in nonlinear Ginzburg-Landau 
equations and quadratically nonlinear optical media, see ^ Q . 

The cubic-quintic nonlinear Schrodinger equation can also been considered a nonrelativistic ver- 
sion of a simple model in nuclear physics that admits nontopological solitons called "Q-balls" 
The results of this paper may be relevant to the issue of the stability of closed cylindrical Q-walls 
raised in ||. 

For radially symmetric waves (to = 0), an energy-based criterion for instability was established 
by Shatah and Strauss |^8| and a stability criterion was established by Grillakis, Shatah and Strauss 
||13| . The result is that waves in a family parametrized by u> are stable if dN / du > and are unstable 
if dN/doj < 0, where TV = L 2 |u 1 2 ■ Consistent with previous studies j|, [l0| |||] we find that all 
spinless waves for the cubic-quintic (1.2) are stable, but if the sign of the quintic term is reversed 
then all waves become unstable. Our numerical results elucidate the mechanism of the transition to 
instability as the coefficient of the quintic term changes sign continuously. A single pair of imaginary 
eigenvalues collides at the origin and a pair of real eigenvalues emerge with opposite sign. A related 
scenario occurs for solitary waves of generalized Korteweg-de Vries equations [ p6| in which unstable 
eigenvalues emerge at the origin (out of the continuous spectrum) as the sign of dN/duj changes, 
and suggests that stability criteria based on the sign of dN/dcu will not detect instability transitions 
that are not associated with eigenvalues moving through the origin. 

For the spinning solitary waves of the nonlinear Schrodinger equation under consideration (m ^ 
0), it is not known whether spectral stability is mathematically a sufficient condition for nonlinear 
stability (modulo symmetries of the wave family as appropriate) . It is conceivable that instability 
could be created by nonlinear resonance phenomena involving both localized oscillatory modes and 
radiation modes in the continuous spectrum, such as have been studied by Soffer and Weinstein 
p9[ . If these waves ultimately turn out to be unstable, however, the absence of unstable eigenval- 
ues suggests that they would enjoy a long "lifetime" under perturbation, and that any instability 
mechanism would be a subtle one. 

This paper is organized as follows. In section 2, we survey the structure of solitary wave solutions 
to (1.1) as th e st anding-wave frequency, spin, and node number vary, contrasting the cubic-quintic 
nonlincarity (|l^) with a pure cubic. We linearize about a solitary wave and prove some fundamental 
results concerning the spectrum of the operator A in section 3. In particular, in section 3.2 we 
establish that we can locate all unstable eigenvalues of A by analyzing a finite number of ordinary 
differential equations. In section 4 we introduce and analyze the Evans functions and Wronskians 
associated with the eigenvalue problem. A representation of the Evans function in terms of exterior 
products, described by Alexander and Sachs j|, proves to be convenient for purposes both theoretical 
and numerical. We present numerical results in section 5. In section 6 we discuss heuristic reasons 
why encapsulated-vortex solutions may become stable in a limiting parameter regime. 



2 Taxonomy of solitary waves with spin 

In this section we illustrate the structure of the solitary- wave profiles w(r) for various choices of spin 
to, node number n and standing wave frequency uj, contrasting two different types of nonlinearity 
for which existence proofs of infinite families of localized standing waves have been given p4| nGl . 



In particular we contrast the cubic-quintic nonlinearity in (1.2) with the pure cubic 

g(u) = \u\ 2 u. (2.1) 



5 



Because of the invariance of (1.1) under spatial reflection, we may assume that m > 0. 



Localized solutions of equation (1.4) for standing wave profiles have the asymptotic behavior (see 



d r m as r — > + , 

w(r) ~ i d x e-° r (2.2) 
' = — as r — > oo, 



for some constants do and doc, where a — yj— f'(0) — y/uj > 0. More precisely, where its amplitude 
is small the behavior of w(r) is governed by the linear equation 



1 



w + -vf T w - a 2 w = 0. (2.3) 

Linearly independent solutions of this equation are the modified Bessel functions I m (ar) and K m (ar), 
which have the asymptotic behavior 



^M~-^->r) m , K m (ar)~2 m - 1 (m-iy.(<jr)- m , (r -» 0+) 



I m (ar) ~ -=L= e°\ K m (ar) ~ e~° r (r - oo) 



2 m m\ 
1 

\/2nar " ' " ' ' V 2or 
except in the case to = 0, when ^(c 7 ") ~ ~ hir as r — > 



(2.4) 



We compute solutions to equation (1.4) numerically by a shooting method, based on small- 
amplitude approximations of the form w(r) ~ doI m (ar) for small r (m > 0), w(r) ~ dooK^ar) 
for large r. For a suitably small initial radius r$, we numerically solve the initial- value problem 



for (1.4) with initial data w(ro) — doI m (aro) and w'(ro) = doaI' m (aro). We vary the parameter 
do (equivalently do) so that at some large radius ri, the solution and its derivative are small and 
assume the proper asymptotic form, satisfying w' (r\) / w{r\) = oK' m {<jr\) j K m {or\). 

In Fig. H we plot the potential F for the cubic and cubic-quintic nonlinearities, taking u> = 1 and 



u) = .175 respectively. We may interpret (1.4) as an equation of motion for a point with position w(r) 
at time r moving in the potential well F(w) — ^m 2 w 2 /r 2 , subject to the time-dependent damping 
force —w'/r. As r — > oo, the particle "balances" on top of the hill at w = 0. Nodes in the radial 
profile occur when the particle crosses the top of this hill. Since the hilltop at the origin is an 
unstable equilibrium of this one-dimensional "motion," solutions satisfying w(r) — > as r — > oo are 
found for do in a discrete set. 

Cubic. In the case of a cubic nonlinearity g{u) — 7|u| 2 u with 7 > 0, we may always scale so that 
7=1 and u> = 1, by replacing w(r) by y^uj~fy w(y/ujr). Figure ^| shows the nodeless (n = 0) wave 
profiles for spin values m — 0, 1, 3 and 6. 

The form of the nodeless waves varies with m in a way that can be understood as follows. 
It appears that as m becomes large, w(r) remains small initially on a large interval, and then is 
approximated near its maximum at r = R by <f> a {r — R), where 4> a is the even positive localized 
solution on (—00,00) of 

(f>" - a 2 cj) + cf> 3 = 0. (2.5) 

with a = y/l + (m/R) 2 . Explicitly, 4> a (x) = \[2 a sech 2 ax. We will argue in appendix A that 
R \/2m, so that max<j) a (x) — v3- This approximation agrees well with the numerical solution 
shown in Fig. |^. 

Cubic-quintic. In the case of a cubic-quintic nonlinearity g(u) = 7|w| 2 it — ^|m| 4 m with positive 
constants 7 and 5, we can scale so there is one free parameter u > 0, and j = 5 = I. The condition 
that F(w) > for some w is then equivalent to the requirement that 

lo < to, = ^ = 0.1875. (2.6) 
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When this holds, F(w) achieves its maximum value at a point w — w max > which depends on u>. 
Potentials F with this additional local maximum are called hilltop potentials in and We 
write F max = F(w max ) and note that u^ ax = \ + — 4uj. For lo = lo* we get w max = a* = v3/2. 
Figure || shows nodeless profiles with lo = 0.18 for spin values m = 0, 1, 2, 3. 

We demonstrate the effect of varying lo in Figures || and ||, where we plot nodeless profiles for 
spins m = 1 and 3, respectively. As lo approaches lo* from below, the height F max of the rightmost 
hilltop in the graph of F decreases to zero, and the "particle" at position w(r) spends more "time" 
near u> max , resulting in flatter profiles. This is the phenomenon of "loitering at the hilltop" studied 
in |H| and |ll| , also see || . It is precisely in this "flattop" regime that we find nodeless waves with 
no unstable eigenvalues for any m, as discussed in Section [|. 

The wave profiles develop a structure in this regime that can be described more precisely. It 



is shown in [[L5| that for any lo G (0,u;*), with an appropriate amplitude d in (2.2) one obtains a 



nonlocalizcd sol ution of (1.4 ) that increases monotonically to u> max as r — -> oo. This yields a vortex 
solution to (1.1). Note that for lo near lo*, the profiles in Figures || and || resemble a vortex profile 
cut off by a kink. Near the kink location at r = R (where w(r) last equals ^w max , say), the profile 
is approximated by (f>*(r — R), where 0* is the solution on (— oo, oo) of 



that satisfies 



0" +.9(0) -u*w = (2.7) 
0*(r) — ► a* as r — ► — oo, <j>*{r) — > as r — * oo, ^*(0) = -a*. (2.8) 



We note that a one-dimensional kink solution of (1.1) of the form u(x,y,t) = e l "**0*( a; ) exists for 



any nonlinearity with h € C 1 that has the properties (i) and (ii) stated in the introduction. 

Fig. ^ suggests also that as the spin m increases, the wave structure contains a flipped inner kink 
that also moves out away from the origin. We formally analyze the asymptotic behavior of the kink 
locations as lo — > ld* in appendix A. 



3 Linearization and reduction to ordinary differential equa- 
tions 



In this section we linearize (1.1) about a given exponentially localized solution of the form (1.3) and 
prove a number of fundamental results that concern the associated spectrum and eigenvalues. In 
particular, we prove in Theorem [TJ] that the search for unstable eigenvalues reduces to the study 
of a finite number of eigenvalue problems for 2x2 systems of second-order ordinary differential 
equations. 

We consider a solution u(x,t) to (1.1) that results from a perturbation to the initial condition 



uq(x,Q) = vq(x) for a standing wave solution uq(x, t) = e vq(x) where Vq\x) 
w(r) — > exponentially as r — > oo. We write 



u(x, t) 



(v (x) + v(x,t)) 



w{r) with 



(3.1) 



and linearize the nonlinear terms in the evolution equation for the perturbation v around amplitude 
v = 0. We compute that 

g(v +v)- g(v ) = a{w)v + P(w)e 2lme v + 0(\v\ 2 ) (3.2) 

where 



"00 



1 



g'{s) + -g(s) 



/?(*) 



1 



s'00 - -g(s) 



(3.3) 



for s / 0, with a(0) — (3(0) — 0. Here g' refers to the ordinary derivative of the real-valued 
restriction of g to real arguments. Note that because g is not holomorphic, its linearization is not 



complex-linear but only real-linear. Also note that g'(s) — 0(s 2 ), so a(s) — 0(s 2 ) and (3(s) — 0(s 2 ) 
as s — ► 0. 

Thus the linearized evolution equation for the perturbation v may be written 

- i d t v = Av - a 2 v + a(w)v + (3(w)e 2im9 v (3.4) 

where a 2 — u> > 0. We next represent the linearized evolution equation as a pair of real equations 
for the components of 

so that v = $i + i$2- Setting 



$2 / \ Im w 



"o 1 )' R =(l -lj> ^ 



we rewrite equation (3.4) as 

9 4 $ = A$ = J (A - a 2 + o(io(r)) + (3(w(r))e 2rn8J R) $. (3.7) 

3.1 Spectrum and eigenvalues 

As a first step toward understanding the stability properties of the standing waves, we investigate 
the spectrum of A, regarded as an unbounded operator on the space L 2 (R 2 , R 2 ) of square-integrable 
functions with values in R 2 . The domain of A is the Sobolev space .ff 2 (R 2 ,R 2 ) of functions into 
R 2 whose derivatives up to order 2 are square-integrable. This space is complexified to study the 
spectrum of A, which we denote by u(A). The spectrum of A consists of two parts, the discrete 
spectrum cr<iisc(^l), consisting of isolated eigenvalues of finite multiplicity, and the essential spectrum 
cr ess (^4), which consists of everything else and is determined by the "operator at infinity" defined by 

A = J(A-a 2 ). 



Lemma 3.1 The spectrum of A is the union of the discrete and essential spectrum. The latter is 
purely imaginary and is given by 

a css (A) = a(A ) = {it \ t e R and \t\ > a 2 }. 



Proof The spectrum of Ag is easily determined using the Fourier transform. It consists of the set of 
all eigenvalues of the matrix J(~£ 2 — a 2 ) as £ varies over R. This set is equal to {it \ t £ R and \t\ > 

The operator A — Aq is a multiplier with coefficients that are continuous and bounded and 
decay exponentially as r — > oo. Consequently, it is easy to show (using the convenient compactness 
criterion of |p5[ , for example) that A is a relatively compact perturbation of Aq, meaning that for 
A ^ o-(Aq) the operator (A — Ao)^ 1 (A — Aq) is compact. From Weyl's theorem, it follows that A and 
Aq have the same essential spectrum. □ 

Remark 3.2 We briefly digress to consider the linearization of ( [L~l| ) about the spatially constant 
solution uo(t) = ae lu)t . One may assume a > without loss of generality, and one has u> = g(a)/a 
and a(a) — u> = (3(a) = a 2 h'(a 2 ). Linearization yields the constant-coefficient operator A a = J(A + 
/3(a) + P(a)R), whose spectrum is the set {±A(£) | £ € R 2 } where A(£) = iy/(\£\ 2 - /3(a)) 2 -/3(a) 2 , 
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which is purely imaginary if and only if (3(a) < 0. Hence the solution ae lujt is spectrally stable if 
and only if h'{\a\ 2 ) < 0. 



By Lemma 3.1, our primary concern becomes the discrete spectrum of A in (3.7). We will search 
for eigenfunctions \E r of A taking values in C 2 . The operator A has real coefficients and has the special 
form A = JL where L is self-adjoint. Consequently eigenvalues occur with Hamiltonian symmetry: 
If A is an eigenvalue with eigenfunction '5 , then A and —A are eigenvalues. This is because '5 is also 
an eigenfunction of A with eigenvalue A, and Jvl/ is an eigenfunction of the adjoint A* = —LJ with 
eigenvalue —A. Note that if A is complex with nonzero imaginary part, then the real and imaginary 
parts of e xt ^> are real solutions to the differential equation (3.7) that are linearly independent. 

The equation A^f = A 'J is equivalent to 



(AJ - a 2 + A + a)tf + (3e 2m9J m = 0. 



Suppose ^ is a square-integrable solution of (3^8). We represent ^>(r,9) € 



1/2 
-i/2 



1/2 

i/2 



(3.8) 

I 2 via eigenvectors of J: 
(3.9) 



where ip+(r,6) and ip-(r,9) are complex- valued. Since 



1/2 

-i/2 



iip- 



1/2 
i/2 



1/2 

-i/2 



1/2 

i/2 



taking the components of 3.8 gives 



(iX - a 2 + A + a)ip + + /3e 2 ""V_ = 
(-iX - a 2 + A + a)i>- + f3e~ 2m0 'ip + = 



0. 
0. 



(3.10) 



(3.11) 



Since and ip- are square-integrable and a(w(r)), f3(w(r)) are bounded, it follows that ip- 
lie in the Sobolev space H 2 (M. 2 ). In particular, they are Holder-continuous functions. 
We next expand ip+ and ip_ in Fourier series with respect to 9: 



(3.12) 



with 



This is an orthogonal expansion and each term lies in H 2 ( 
to shift the summation index and write instead 



m0 d9. (3.13) 
2 ). It will prove notationally convenient 



J(j±m)9yU±m) 



V). 



(3.14) 



The Fourier coefficients then satisfy 

{iX -a 2 + A r - r- 2 (j + mf + a)^ +m) 
(-iX -a 2 + A r - r- 2 (j - m) 2 + a)y U ~ m) 

and satisfy the boundary conditions 

lim y± ±m \r) exists, and lim y± 

r — >0+ r — 



{j±m) 



(r) = 0. 



0, 
0, 



(3.15) 



(3.16) 



Here A r = 4-$ + - jj^ is the radial Laplacian. 



- ' r" 1 

For each integer value of j, equations (p. 15) are a system of second-order ordinary differential 



equations for the two functions y^_ +m ^ (r) and y K _l m> (r). There is no coupling between these systems 
for different values of j. Summarizing, we have proved the "only if" part of the following: 

Proposition 3.3 A complex number A </ a css (A) is an eigenvalue of A in L 2 if and only if for some 
integer j , equations (3.11 ) have a nontrivial solution satisfying (3.16). 



T, 



To prove the "if" part, we claim that if there are nontrivial solutions y± to equations 



( |3.15 )-( 3.16 ) for some A and some j, then 



1/2 
-i/2 



(r) 



1/2 

i/2 



(3.17) 



is an eigenvector of A with eigenvalue A. Clearly ^ = \I/<J' A ) satisfies (3.8) except perhaps at the 
origin. We must show that iff S H 2 (R 2 , M 2 ). It follows from (2.2) and the theory of asymptotic 
behavior of ODEs (see and Subsection 4.1 below) that solutions of ( 3.1 5| ) that decay to zero as 
r — > oo must decay at an exponential rate, together with their derivatives. Therefore iff S L 2 (M. 2 , M 2 ). 
Furthermore, solutions that approach a limit as r — > must have bounded derivatives. It follows 
that iff satisfies (3.8) in the sense of distributions, and from this it follows that iff € -ff 2 (R 2 ,R 2 ). 



This proves the Proposition. D 

Corresponding to the complex eigenvector (3.17) of the operator A are the two real- valued solu- 
tions 

$£ A) (<) = Re(e At ^'' A) ), $[ J fc ) A) (i) = Im(e At ^'' A) ) (3.18) 



of the vector evolution equations (3.7), from which we recover two complex- valued solutions of the 
original linearized equation (3.4) via 



«W(t) = *g$ ) (t)+»^(t) 
Explicitly, we have: 

v (a *>(t) = I (e xt e i ^ +m ^ e y + + e Xt e i( j-™)<>y_^ , 

Ti (' 



(*) = ^(*) + i^(*)- 



e \t e i(j+m)O y+ 



e \t e i(j-m)0y_ 



(3.19) 



(3.20) 



Since any real linea r co mbination of vr a -\ and vn,\ is also a solution, the general solution to the 
linearized equation (3^4) corresponding to eigenvalue A and index j is 



v(x, t) = (ce xt e^ +m ^y { l +m ' X) (r) + ce ^U-m)e y U-m,X) ^ ^ ^ ^ 

where c is any complex constant. If wc were to take the time-zero value of this function v(x, 0) as a 



perturbation to the initial conditions Vq{x) for (1.1), we would expect the nonlinear time evolution 
of v to be well-approximated by ( |3.21[ ) for small time. 

We note that j measures the departure of the perturbation v from the angular dependence of 
uq. We will refer to j as the "twist" index. 



3.2 Bounds on index for unstable eigenvalues 

We have shown eigenvalues of the operator A correspond to eigenvalues of ( 3.15 )-( |3.16|) for some 
value of the index j. To locate eigenvalues of the operator A, we use numerical techniques to locate 
eigenvalues of ( 3.15 )— ( 3.16 ). This procedure at first appears intractable because there is a countably 
infinite set of equations to analyze. Our main result in this section, however, establishes that any 
eigenvalue with nonzero real part occurs in a system (3.15) with twist j in a finite set determined 
by the wave profile w alone. 
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Theorem 3.4 Suppose that Re A ^ 0, where A is an eigenvalue of A corresponding to a solution of 



(3.11 )—(3.16) for some integer j. Then 

\j\ < imax := \m\ + \J~C~i, where C 2 = maxr 2 (a(w(r)) + \j3(w(r))\) . 

r>0 

Proof Let fi(s) = a(s) + \(3{s)\. Then ^(0) = and fj,(s) = cr 2 + max{/'(s),/(s)/s} for s ^ 0. 
Because / is continuously differentiable, /i is continuous on R, and we note that ^t(r) takes on positive 
values because f(w(r)) does (otherwise the primitive F(w(r)) would never get positive). Because 
n(s) = O(s) as s — * and w(r) is exponentially decaying at oo, the maximum C 2 of r 2 /i(u>(r)) is 
achieved, and we have < C 2 < 00. 

We note that with y± = y£ ±m \ (3.15) can be written 



where 



and 



L jY = XKy, (3.22) 



_ ,'a 2 -A r + r- 2 (j + m) 2 -a -(3 \ , . 

Lj ~\ -13 a 2 - A r + r- 2 (j — m) 2 — a J ' {6M) 



Define inner products in L 2 (M. + ; r dr) and L 2 (K + , R 2 ; r dr) by 

/>oo 

(y,z)= y{r)z{r)rdr, (y|z) = (y + , z+) + z_). 



We claim that if C* := min{(j — m) 2 , (j + m) 2 } > C 2 a nd (3.22 ) has a nontrivial solution, then 
(Ljy\y) > 0. To see this, recall that for a solution of ( 3.15 )- fl3.16| ), y± and its derivative are bounded 



as r — > and decay exponentially as r — > 00. We then have that ((cr 2 — A r )y_|-, y±) > 0, therefore 



DC: 



U iy \y) ■- I U J — \ 2 \y + \ 2 +( J —^) 2 \y-\ 2 )rdr 



r 



> 



3C 



{" (I2/+I 2 + \y-\' 2 ) + Hv+y- +y+y-)}rdr 

(C,-C 2 )(\y + \ 2 + \y^\ 2 )r- 1 dr>0, 



where we used the facts that yz + yz < \y\ 2 + \z\ 2 and a + \/3\ < C 2 /r 2 . 

Now take the inner product of ( 3.22| ) with y. Since C* = (|j| — |m|) 2 , if \j\ > \m\ + \fC2~ we get 



0< (L,y| y ) =\(Ky\y). 
Since K is skew-adjoint, (Ky\y) is purely imaginary, therefore A must be purely imaginary. D 

3.3 Bounds on eigenvalues 

We show here that eigenvalues of A are confined to a vertical strip determined by the wave amplitude 
profile w. 

Proposition 3.5 If X is an eigenvalue of A, then 

I Re A| < A max := max(\a(w(r))\ + \p(w(r))\) . 

r>0 V / 
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Proof Let A = A-A Q = J(a + f3e 2mej R). Since the I 2 matrix norms \J\, \R\ and \ e 2mf>J \ all equal 
1, we have \A\ < Cq := max r .>o(|a;(w(r))| + \/3(w(r))\). Note that A — Ag is the Fourier multiplier 
F(M) defined via (.F(M)#)"(£) = M (£)*(£), where 



M(0 = 



2 _ U |_> 

A 



-«7- - |f| 



(3.25) 



Suppose $ is an eigenfunction of A corresponding to an eigenvalue A with nonzero real part. 
Then = (A - A - A)9, and A - A = T{M) is invertible, so that we have 



(7 — J r (M~ 1 )A) i $> = 0. 
Let V± = A ± i(a 2 + |e| 2 ). Then A = \{V+ + V-), a 2 + |£| 2 = ±{V+ - V-), 



(3.26) 



M(£)=QdiagfP-,7>+}Q~\ Q = 



= J_fl 1 



Note that Q is unitary and M^y 1 = Q diag{"P_ 1 , V + 1 }Q~ 1 , so the matrix norm |M(^) _1 | < 
maxlT 5 ^ 1 ! < 1/| ReA|. Therefore, the operator norm o f J 7 (M ~ 1 )A on L? is bounded by Cq/\ Re A|. 
If | Re A| > Co, then, there is no nontrivial solution of ( 3.26 ) in L 2 . □ 



3.4 Symmetries 

The value A = is always an eigenvalue of the operator A, with multiplicity greater than one, due to 
various symmetries of the original nonlinear Schrodinger equation (e.g., translation, rotation, 

phase change). These symmetries and associated generalized eigenfunctions are described in detail 
in Appendix M. 



The eigenvalue equations (3.15) also have simple symmetries connected to the Hamiltonian sym- 



metry of A. From (3.15) and the orthogonal expansion in (3.14) it is easy to check the following 



Lemma 3.6 Suppose form a solution to ( 3.15 ) for some pair Then (y-,y+) form 

a solution for (— A, — j) and (y_,y + ) form a solution for (A, — j). If X is a real eigenvalue of A 
corresponding to index j ^ 0, then X is (at least) a double eigenvalue. 



4 Detecting eigenvalues with Evans functions 

In this section we prove that eigenvalues of the operator A of index j correspond to zeros of a certain 



analytic function Ej(X), which we call an Evans function for the system (3.15). This characterization 
of eigenvalues as zeros of an analytic function allows us to use the argument principle to count the 
total number of eigenvalues inside various contours in the complex plane. 

We shall provide three equivalent expressions for this Evans function, in terms of a Wronskian, 
an adjoint system, and an associated exterior system. The last expression is used to establish the 
global analyticity of Ej(X) and is advantageous for reasons of numerical stability. 

We develop the theory in this section only for the case of spinning waves for which m/0. The 
spinless case m = is different in several ways, primarily due to the nonzero value of the wave profile 
w(r) at r — 0. For brevity's sake we omit any theoretical treatment of the spinless case, though in 
section 5 we report results of numerical computations performed for this case. 

4.1 Asymptotic behavior of eigenfunctions. 



Here we study the asymptotic behavior of solutions to ( 3.15 ) as r — > + and as r — > oo. We fix the 
integer m ^ 0, the base wave radial profile w, the index j, and a complex number A in the cut plane 
C a = C \ {ir : \t\ > a 2 }. 
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Because w(r) — > exponentially as r — * oo, for large values of r the coefficients a(w(r)) and 
/3(w(r)) in ( 3.1 5| ) vanish exponentially quickly as r — + oo. We thus expect that solutions to ( 3.1 5| ) 
will behave for large r like solutions of the decoupled equations 



(£A - a 2 + A r - r~ 2 (J + m) 2 )y + = 0, 

2 



(4.1) 

(-iA -a 2 + A r - r~ 2 (j - m) 2 )y_ = 0. (4.2) 

Solutions of ( |4.lD are linear combinations of the modified Bessel functions I\j+ m \ (k+r) and K\j+ m \ (k+r), 
and solutions of ( |4.2j ) are linear combinations of the modified Bessel functions I\j_ m \(k_r) and 
K\j_ m \(k-r), where 

k 4 



(4.3) 
(4.4) 



\/ a 2 — iX, k- = \/ a 2 + iX, 
with the principal branch of the square root taken. For brevity we write 

7 +( r ) = I \]+ m \{k+r)-, I-(r) = I\j- m \(k-r), 
K+(r) = K lj+ml (k + r), K^{r) = K\j_ m \{k-r). 

Since A is in the cut plane, the complex numbers k± have positive real parts, and thus |i±(r)| — > oo 
and |iT±(r)| — * as r — > oo. 

Similarly, because w(r) — > as r — > + when m > 1, we ex pec t that solutions to (3.15) will also 
behave in this case like solutions of the decoupled equations ( |4.l[) , (1.2) for r near zero. Note that 
|l±(r)| is bounded and |2f±(r)| — > oo as r — * 0. 



To make these observations precise, we reformulate ( 3.15 ) as a first-order system. We set 



c+(r) = fci + 



j + m 



c_(r) = fc 2 + 



3 — m 



Then (3.15) is equivalent to the system 

y' = B(r,i, A)y 
with B(r,j, X) — B(r,j, X) + B(r), where, writing y + = y+ r ' 



= V- 



, we have 



y 



(V+\ 




vA 


'-( 


V- 




\V-'J 





1 

-1/r 









c_ 



B = 



0^ 

-a -P 



-P -a 0, 



(4.5) 



(4.6) 



(4.7) 



Let y i5 i = 1, . . . , 4, denote the linearly independent solutions of the asymptotic system 

o 

y' = B(r,j, X)y 
that form the columns of the fundamental matrix 

'J+(r) if + (r) 

j;(r) o jr;(r) 

J_(r) 
I'_(r) 




(4.8) 



(4.9) 



Lemma 4.1 fix A m the cut plane C a — C \ {it : |r| > a 2 }, and let m > 0. Then there exist 
solutions and y\°°^ to the coupled system such that 



y\°\r) 



.(oo) 



(r) 



y { (r) as r 
yi(r) as r 



0+ 
oo, 



(4.10) 
(4.11) 



fori = 1,2,3,4. 



The proof of Lemma 4.1 is in Appendix C.l 
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4.2 Evans function test for eigenvalues 

We fix the spin m > 1, base wave w, index j, and a complex number A in the cut plane C CT . To 
determine whether A is an eigenvalue corresponding to index j, we wish to determine whether ( 3.1 5| ) 
has a nontrivial solution satisfying ( 3.16 ). This is equivalent to determining whether ([0]) has a 
bounded solution that vanishes as r — > oo. 

1, 

From Lemma 



,4, be solutions of (4.6) with the pro perties specified in Lemma 4.1. 
,4.1 we know that the subspace So of solutions y to (4.6) that are bounded as r — * + is 
spanned by y fcind y^ ■ Also we know that the subspace S, 
oo is spanned by and y± 



of solutions y to (4.6) that approach 

zero as r — * 00 is spanned by y 3 °°^ and y^°\ 

The question of whether there are nontrivial solutions of ( 3.1 5| ) satisfying (3.16) for the value 
of A under consideration is thus reduced to the question of whether the subspaces So and Sx, have 
nontrivial intersection. To answer this question, we may compute the Wronskian of the four solutions 
yf \ y^\ 1/3 an d 2/£ • Since there is a nontrivial intersection of So and Soo if and only if these 
four vectors are linearly dependent, the Wronskian of these four vectors vanishes (identically in r) 
if and only if the number A is an eigenvalue corresponding to index j. 
We note that, since tr(B) — — 2/r, by Abel's formula the quantity 



Met 



(00) 
2/3 • 



(00) 



(4.12) 



is independent of r. This normalized Wronskian is called the Evans function [El associated with the 



system (4.6). Thus we have the following. 



Proposition 4.2 A complex number A in the cut plane C CT is an eigenvalue of A if and only if 
Ej(X) = for some integer j. 

We conclude this subsection by describing some symmetries of this Evans function. 

Proposition 4.3 For all integers j and complex numbers A 6 C CT , 

(i) E^Xj = E j (-\) 

(ii) Ej(X) is real if X is purely imaginary, 
(tit) Ej(X) = E-j(—\) 



Proof (i) Since taking the complex conjugate of (4.6) yields a solution of the same equation with 
A replaced by —A, we find from (4.12) that Ej(X) — Ej(—X). 
(ii) follows from (i), since A = —A if A is purely imaginary. 



(iii) We mak e use of the symmetry of equations (3.15) as described in subsection 3.4. If y is a 
solution of (|4.6| ), then exchanging the first two components with the last two yields a solution of 
y' = B(r, — j, —X)y. Consequently, E_j(— X) is obtained from Ej(X) in ( 4.12 ) by exchanging the 
first two rows with the last two. This leaves the determinant invariant, and then (iii) follows. D 



4.3 The adjoint system 

To obtain alternative characterizations of the Evans function, we consider next the adjoint system 
associated with (4.6), 

z' = -zB(r,j,X). (4.13) 

Here z(r) is considered to be a row vector. Any solutions of ( 4.1 3| ) and (4.6) have the property that 
the product z(r) ■ y(r) is independent of r. The asymptotic adjoint system 



-zB(r,j,X) 



(4.14) 
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has the fundamental set of solutions z^, k = 1, . . . , 4 given by the rows of the matrix 



/ n n _ V !<v\ 

Z(r) 



r + (r) -I+(r) 








(4.15) 



I'_(r) 

o o o 

This matrix is the inverse of Y(r), i.e., we have ZY = I. 

Lemma 4.4 The adjoint system has solutions z^°\ k — 1, ... ,4 that satisfy 

(oo)/ n / \ 

z k '(r) ~ Zk(r) as v — ► oo . 



This lemma is proved in the same manner as Lemma 4.1, so we omit the details. 

Let z k °° , k — 1, ... ,4, be solutions with the properties specified in Lemma i.4. Let Z^-°°\r) 
denote the matrix whose rows are z k °°\r), and let Y^°°\r) be the matrix with columns y k °°\r). 
Because of the limiting behavior as r — > oo expressed in Lemmas 4.1 and 4.4 and the fact that 

o o 

ZY = I, it follows directly that 

2H( r )yW( r ) = j (4.16) 

for all r > 0; that is, z { k °°\r) ■ y, (oo) (r) = S kl . 

The matrix y(°°) is a fundamental matrix for the system (4.6). Therefore the solutions yf^ of 
Lemma 4.1 admit representations of the form 



(o) 

with constant coefficients that are given by 

Oiki 



E(oo) 



(4.17) 



fe=i 



The value A is an eigenvalue of fl3.15| )-(3.16) if and only if the span of y± and y 2 °^ nontrivially 



.(°) 



intersects the span of y^"^ and y^' ■ This is equivalent to the statement that 



det 



( an 


au ^ 


= det ( 


V "21 


"22 > 





(°°) ..(o) „(«>) ro) 



■ - Z 2 ' 2/l Z 2 

It turns out that this determinant simply yields an alternative characterization of the Evans function. 
Proposition 4.5 For all integers j and complex numbers A S C CT , 



(0) 



„(«) 

'2 



y 2 



o. 



(4.18) 



£j-(A) = det 



(oo) „,(0) (oo) (0) 



• y 



•1/2 



,(°°) „,(«) „(») „,(°) 



Proof 



3/i 

,(°o) .,(<») 



• V2 



= -r*det (y^U^U' 

4 

= -r 2 £ aiia fe2 det ( y| M \ y£ , y|™\ yi ^ 

i,k=l v 

9 / \ i , / (oo) foo) foo) (oo) 

= -r z (a u a 2 2 -a 2 iai 2 )det ly^ ; ,y 2 ',y^ %y£ 

= -r 2 det fy^°°^ (ana 2 2 - a2i«i2) 
= -r 2 det (V (0) ^ (a n a 22 - a 2 i«i 2 ) 



det 



(oo) (0) (oo) (0) 



■ y 



■y 2 



(oo) „,(0) Joe) (0) 



y 2 



Here we have used the fact that the r-independent quantities r 2 det \ V^ 00 '^ and r 2 det ^Y"( ^ 
are equal by virtue of (4.11). D 
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= -1 



4.4 Exterior systems and analyticity 

According to our results so far, for each fixed A in C CT , Ej(X) depends upon a choice of certain 



.(0) 



and y^°°' > or zf" 1 having the properties described in Lemmas [4.l| and |4.4|. In this 



solutions y 

section, we establish that the value of Ej(X) is independent of this choice, and that Ej is an analytic 
function in the cut plane <C a . 

The method of proof of Lemmas 4.1 and 4.4 could be used to establish that the required solutions 
can be chosen to depend analytically on A, locally in the set C CT in the complex plane. However, 
this approach does not guarantee that there is a single globally analytic choice, leaving the difficulty 
that a given set of choices may yield an Evans function that is not analytic in all of C a . 

This difficulty is obviated with an algebraic device described by Alexander and Sachs || (used by 
one of the present authors for numerical calculations mentioned in that paper) and recently employed 
by Afendikov and Bridges 0. A similar device (called the compound matrix shooting method) has 
been used in numerical computations for boundary- value problems related to hydrodynamic stability 
theory pi. 



One considers exterior products of solutions of the system (4.6) and the associated adjoint system 



(4.13). In the present context it suffices to consider exterior products of two solutions at a time. 



In general, suppose y and y are vectors in C™ that, with respect to the standard basis vectors 



ej, j — 1, . . . , n, have the expansions 



y = ^yj e i> y 



E 



(4.19) 



Let N = 



l)/2 and let y A y denote the vector in C with components 



(V A y)j A k = VjVk - VkVj = det ( V j V J 

\ yk yk 



(4.20) 



for j < k. The vector y Ay = — y A y is called the exterior product of y and y. For convenience, 
when n = 4we order the components according to the identification 

(1A2, 1A3, 1A4, 2A3, 2 A 4, 3 A 4) = (1, 2, 3, 4, 5, 6). 

Suppose that y and y are both solutions of the linear system 

y' = By (4.21) 

where the matrix- valued B has columns bj, j — 1, . . . ,n. Then since By = ^2l/jbj, one computes 
that y := y Ay satisfies the linear system 



y 



(B A )y 



(4.22) 



where matrix B A has columns given by bj A k — bj A + ej A bk- The system in (4.22) is called the 
exterior system associated with (4.21). Similarly, for two solutions z and z of the adjoint system 



z' = -zB 

the exterior product z := z A z is a solution to the associated exterior system 



z' = -z(B A ). 



(4.23) 
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It is a straightforward matter to compute that 



[z A z) ■ (y A y) = det 



z ■ y zy 

zy zy 



(4.24) 



As a consequence of (4.24) and Proposition 4.5, the Evans function Ej may be expressed as 



E j (\) = (z^Azt ) )-(y? ) AV$' ) ). 



(4.25) 



This expression has advantages for both theoretical and numerical purposes, which stem from the 



fact that the exterior products A z£"' and y^' A y^' have a separate characterization that 
shows they are independent of any particular choice of the individual solutions involved. 

Proposition 4.6 Fix an integer m ^ 0, the base wave radial profile w, the index j , and a complex 
number A £ C a . The exterior system associated with (4-6) has a unique solution y^°\r, X)with the 
asymptotic behavior 

y K ' ~ 2/i A y 2 as r -> 0+. 

Similarly, the exterior system associated with the adjoint system (4-1S) has a unique solution 2^°°'(r, A) 
with the asymptotic behavior 

-Coo) ° ° 

z ' ~ z\ A Z2 as r — ► oo. 
Furthermore, the solutions y {r,X), z^°°'(r, A) are analytic functions of X globally in the cut plane 



v (°°) 



(0) 



(0) 



Cor olla ry 4 .7 For any choice of solutions y^\ y^°\ z k having the properties indicated in Lem- 
mas l.l and (.4 we have 



(0) . (0) ~(0) 

y\ ^y 2 =y , 



z^ Az^ = z (oo) . 



The Evans function Ej is given by 



Ej(X) 



j(°°) . y(0) 



(4.26) 



and Ej is analytic globally in the cut plane C„ 



The proof of Proposition 4.6 is in Appendix C.2. The uniqueness property of the solution y 



,-,(0) 



with the given asymptotic normalization stems from the fact that it is a solution of the exterior 
system associated with ( |4.6| ) having maximal decay rate as r — * 0. Similarly, z^ 00 ^ is a solution of 
the adjoint exterior system having maximal decay rate asr-> 00. 



Corollary |4.7| follows from Proposition [4.6| via the observation that y± Ay^ and z^' Az^' are 
solutions of the exterior systems associated with (4.£) and ( f,13| ), respectively, with the asymptotic 



behaviors given in Proposition |4.6[ Thus both yf^ A y£' and z^ 1 A z^' are globally analytic 



,(0) 



,(00) 



„(°°) 



functions of A in C a , and it follows from ( 4.25| ) that Ej is given by (4.26) and is also globally 
analytic in C CT - 



5 Numerical computation and results 
5.1 Evaluation of Evans functions 



To evaluate the Evans function Ej(X) in (4.12) for given j and A, we use a Runge-Kutta ODE 
package (rksuite) to compute numerical approximations to the solutions y^ir) and yf\r) of (4.6) 



that are bounded at the origin. To do this, we solve two initial value problems for the system (4.6) 
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on an interval r$ < r < r^, where ro and are chosen so that w(r) is small (usually < 10 4 ) on 
(0, ro) and (V^, oo). Then a(w(r)) and /3(iu(r)) are close to their asymptotic values on the intervals 
(0,ro) and (^,00). We make the approximations y^iro) 2/i(ro) (* = 1,2) to determine the 
initial conditions. We then evaluate the Wronskian of the computed approximations to y^^oo) 
and y 2 °' '( r oo) with the explicitly known approximations y^ijoo) and y 4 (roo) to the respective values 
y3°°^( r oo) and j/4°°' > (joo) of the solutions to the system ( p~6| ) that vanish at infinity. 

Our actual numerical implementation of this scheme makes use of a change of variables and scaling 
to reduce the sometimes rapid variation of solutions and improve the quality of the computation. 



We represent a solution y(r) of the system (4J3) as an r-dependent scaled linear combination of the 
basis elements |y, : i = 1, . . . ,4| for the asymptotic system (^8), writing 

y(r) = Y(r) exp(Kr)b(r) (5.1) 

where K :~ diag{— k + , — k + , fc-}, Then we can anticipate that the vector b(r) is well-behaved 
as r becomes large, and it satisfies the system 

b'(r) = D(r)b, (5.2) 

o _ o 

whose coefficient matrix D(r) = K + e~ Kr Y~ 1 (r)B(r)Y(r)e Kr remains bounded as r — * 00. 

The Evans functions Ej(X) are also approximated in a similar but simpler way based on the 
expression in ( 1.26| ), by computing an approximation to e~ kr y^ (r) for r < r < and using 



the known asymptotic approximation for ^"'(foo). The scaling factor k is given by k = k + + fc_. 
Numerical convergence studies indicate that we usually attain 8 to 10 digit accuracy with this 
approach. This method of evaluating Ej(\) is very stable, since y^ is a solution of the exterior 
system with maximal asymptotic growth rate. It also avoids a potential difficulty one can encounter 
with the Wronskian approach — maintaining the linear independence of approximations to and 
$4 i s difficult if the growth rates for these solutions differ greatly and if the interval of computation 
is large. 

5.2 Counting eigenvalues with the argument principle 

Let us first recall what is known about the possible location of unstable eigenvalues. According to 



Proposition 3.3 we may locate eigenvalues of A by locating zeros of the Evans functions Ej(X). For 
a given standing wave uo — e lut e lme w(r), Proposition |3.5| guarantees that all such zeros lie in the 
strip where 

I Re A| < A max := max(|a(«;(r))| + \0(w(r))\). (5.3) 

r>0 



Theorem 3.4 guarantees that all unstable eigenvalues A will be zeros of Ej with 



|j| < imax := H + max ^fr 2 (a{w{r)) + \(3(w{r))\) . (5.4) 

r>0 



Furthermore, it follows from Proposition 4.3 that E-j(X) — if and only if Ej{\) = 0, so that zeros 
of E_j in the strip < Re A < A max are in one-to-one correspondence with zeros of Ej in that strip. 

To search for zeros of the Evans functions in the strip | ReA| < A max , we apply the argument 
principle to count the number of zeros of the analytic functions Ej inside a contour that encloses a 
"large" part of the strip. We do not have a rigorous bound on | Im A| for unstable eigenvalues, but 
in practice we find that \Ej(X)\ grows for large |AJ, so we take a large contour determined so that 
an asymptotic trend appears established. Figure [7] shows a schematic of the contour 7 employed in 
our computations, lying in the cut plane C a - The contour is designed to enclose the region where 
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| Re A| < A max and | Im A| < K, except for an e-neighborhood of the essential spectrum {it \ \t\ > a 2 }. 
The parameter e is taken small (~ 10 -6 ), and K large. 

Using the computed approximation to Ej{X) for given j and A described above, we compute the 
winding number of the image £^(7) about 0. This is the change in the argument of the complex 
number Ej(X) as A traverses the curve 7, divided by 2tt. We use an adaptive stepping procedure 
to traverse 7 in steps that result in a small relative change in Ej(X), and accumulate the change in 
aigEj{\). By the argument principle, this yields the number of zeros of Ej enclosed by 7. 



Because of the symmetry Ej(— A) = Ej{\) from Proposition 4.3 (iii), in practice we need only to 
compute the change in argument of Ej{X) along the portion 7,, of the contour 7 in the right-half 
plane Re A > 0. Doubling yields the total change in argument along 7. 

By the same symmetry, the number of zeros enclosed by 7 in the right half-plane is the same 
as the number in the left half-plane. So, from the winding number above, we need to subtract the 
number of purely imaginary zeros and divide by two to determine the number of unstable eigenvalues 
of A corresponding to index j enclosed by 7. To find the purely imaginary zeros we exploit the fact 
established in Proposition that Ej (A) is real if A is purely imaginary. We plot the real function 
Ej(it) for t £ [~u 2 + e, cr 2 — e], refining and checking the multiplicity of zeros with log- log plots as 
necessary. 

The computation of the winding number can be considerably slowed by the presence of zeros of 
Ej(X) very near the contour — the adaptive stepping procedure takes small steps to track the image 
curve near such a zero. The reason we do not take our contour to be a simple rectangle with one 
side along Re A = e is to avoid passing near numerous zeros on the imaginary axis between — ia 2 
and ia 2 . 

Nevertheless, it sometimes happens that the values of Ej(X) become small as A traverses the 
part of the contour j r in the right half-plane near the essential spectrum on the imaginary axis. 
This behavior suggests that the contour passes near a zero of Ej . To verify that these zeros do not 
lie in the right half-plane outside 7, we use low-degree least-square polynomial fits to extrapolate 
an approximation to Ej(X) near j ri and a root finder to locate the zeros. In all cases the zeros 
appear to lie in the left half-plane. In view of the symmetry Ej(X) — Ej(—X), we interpret these 
findings to indicate not that Ej has zeros in the cut plane C a outside 7, but rather that the analytic 
continuation of E-j(X) across the imaginary axis has zeros near the axis. Such points are known as 
resonance poles, see ]26| . 

5.3 Spinless ground-state waves 



The stability of ground state (nodeless) standing- wave solutions to (1.1) with nonlinearities that 
include those we consider here is analyzed in jl3| and |54|. These works establish that a family of 
ground state standing waves (m = n = 0), parametrized by standing wave frequency ui, is orbitally 
stable (stable modulo spatial translations and phase shifts under H 1 perturbations of initial data) 
if dN/duj > and unstable (to radial perturbations) if dN/duj < 0, where N — ||wo|| 2 - For general 
nonlinearities, the stability of the ground state in the marginal case dN/duj — is open, but for pure- 
power nonlinearities g(u) = r y\u\ p ^ 1 u (for which dN/duj = is equivalent in d spatial dimensions 
to p = 1 + 4/<i) the ground state is unstable in the sense that there exist arbitrarily small L 2 
perturbations of initial data for which the corresponding solution of (NLS) blows up in H 1 norm in 
finite time. (See @.) 

The construction and computation of the Evans functions for these spinless waves differs from 
the description we have given above, since the wave amplitude w(r) approaches a nonzero value as 



r approaches zero. Nevertheless, through diagonalization of the coupling matrix in (3.15), explicit 
(formal) asymptotic solutions may be written and Evans functions computed. For brevity we omit 
the details, since our primary concern is the spectral stability of standing waves with spin. 

We performed computations for spinless waves for the cubic and cubic-quintic nonlinearities, 
both to validate the numerical code and to study the mechanism of transition to instability as the 
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quantity dN/duj changes sign. 

The cubic nonlinear Schrodinger equation in two spatial dimensions is a marginal case from the 
point of view of theory (dN/duo = 0). Recall that by scaling we may without loss of generality 
take g(u) = \u\ 2 u and consider only the fixed standing wave frequency u> = 1. The bounds from 
Theorem [O] and Proposition 3.5 yield that |j| < j max = 1.85 and |ReA| < 14.6 for unstable 
eigenvalues. In fact we find no eigenvalues with positive real part, indicating that the mechanism 
for instability in this case is not exponential linear instability. The eigenvalue A = has high 
multiplicity, however — it is found to be a zero of Ej of order 4 for j = and of order 2 for j = 1. 
(This remains true for the cubic for all spin/node- number combinations considered.) This order is 
consistent with the enumeration of zero modes in Appendix B. No other discrete eigenvalues are 
found. 

To study transition to instability, we perturb the cubic nonlincarity and analyze the one-parameter 
family g(u) — \u\ 2 u — 5|w| 4 w. Making use of the scaling transformation w(r) — y/ujw(y/ujr), we may 
without loss of generality study eigenvalues for the single standing- wave frequency uo = 1. For the 
wave to exist, 5 must satisfy S < j^. Numerical computation of N = 2n J °° w(r) 2 rdr for the wave 
profiles indicates that, by the stability criterion for ground states, uq is stable (dN/duo > 0) for 
< 5 < ^ and is unstable (dN/duo < 0) for 5 < 0. 

Our numerical search for eigenvalues shows that the transition to instability of the ground state 
as 6 decreases through zero is characterized by the collision of two purely imaginary conjugate 
eigenvalues in the twist-0 (j = 0) subspace. We observe the transition by plotting E (ir) for 
r S (—a 2 , a 2 ) for various values of 8. For < 5 < f^, log-log plots indicate that the function E (ir) 
has a zero of multiplicity two at r = and two other conjugate simple zeros. As S — > + , the simple 
zeros approach r = 0, and at 6 = there is a zero of multiplicity four at r = 0. For 6 < 0, two real 
zeros of Eo(X) emerge with opposite sign on the real axis. This creates a positive eigenvalue that 
renders the standing wave exponentially unstable. 



5.4 Waves with nodes 

In all our computations for "first excited states," waves whose profiles have a single node where 
w(r) — for some r > 0, we find many unstable eigenvalues. This holds whether the waves have 
spin or not. For example, for the cubic with m = 0, n = 1, we find 12 unstable eigenvalues, 
distributed as (2,1,1,1,1,1,0...), which we abbreviate as (2, l(5x ), . . .), meaning two in index 
j = and one each in index j for |j| = 1, ... ,5. For m = 1, we find 24 unstable eigenvalues, 
distributed as (2(4x ), l(5x ), . . .). For m = 2, we find 34, distributed as (2(6x), l(6x),0. . .). The 
existence of eigenfunctions with twist j = leads us to expect that the higher-node-number states 
have exponential instabilities to perturbations of initial data of the form e" n6 v(r). For the cubic- 
quintic g(u) = \u\ 2 u — |m| 4 m we only mention a few cases with n = 1 since now different values of 
ui are in principle different: with m = 0, to = 0.05 we find 10 unstable eigenvalues; with m = 1, 
uj = 0.16 we find 20; and with m — 2, oj — 0.165 we find 24. It proved difficult to perform accurate 
computations for profiles with more than one node (n > 1), but we expect all such waves to be 
unstable. 



5.5 Spinning waves with no nodes 

For the cubic nonlinearity (scaled so uj = 1) our numerical results indicate that spinning standing 
waves with no nodes (n = 0) are unstable. With spin m = 1, there are 6 unstable eigenvalues^], one 
for each twist index with |j| = 1, 2, 3. And there are 8 eigenvalues in the "gap" {it \ —a 2 < t < a 2 }, 
distributed in j as (4, 2, 1,0,1,0.. .). For m = 2, there are 6 unstable, one for each j = 1, . . . , 6, 
with 11 gap eigenvalues distributed as (4, 2, 1, 1, 1, 0, 0,1,1,0...). 

l 1.17 ± 1.52i, 1.97 ± 1.69i, 1.81 ± 1.88* 
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Table 1: Eigenvalue distribution and bounds at the instability transition for the cubic-quintic non- 
linearity. Waves with spin m are spectrally stable for u> cr < uj < uj*. Unstable eigenvalues emerge 
at A cr , and can exist only for |j| < j max by Theorem 3.4. 



m 




A C r 


jmax 


=ff of zeros in gap, for j = 0, 1, . . . 


1 


0.1487 


0.0478i 


7.66 


4,3,3,2,1,0... 


2 


0.1619 


0.027U 


14.2 


4,4,3(4x),l(4x),0... 


3 


0.1700 


0.0136i 


23.4 


4,4,5,5,4,3(4x),2,l(5x),0... 


4 


0.1769 


0.0063i 


38.3 


6(5x),5(4x),3(7x),2,l(7x),0... 


5 


0.1806 


0.0033i 


57.6 


8(3x),7(6x),6,5(4x),3(10x),2(3x),l(10x),0.. . 



Next we consider the cubic-quintic nonlinearity g(u) — 7|it| 2 w — <5|it| 4 w, where the constants 7 
and S are both positive. As mentioned in section ||, the scaling w(r) — w(-^=r) allows us to 
suppose without loss of generality that 7 = £ = 1. The allowable range for u>, for which F(s) > 
for some s > 0, is < uj < uj* — j^. 

Spectrally stable waves. For the nodeless waves with spins m — 1,2,3,4,5, we have discovered 
ranges of standing- wave frequency uj for which all the Evans functions Ej, j = 0, 1, . . . , jmax, appear 
to have no zeros with positive real part. Our numerical computations indicate that, for a given 
value of spin m, there is a transition frequency w cr such that A possesses unstable eigenvalues for 
< u> < lu ci but A has no unstable eigenvalues for w cr < uj < uj* . In Table [l] we list the transition 
frequency uj ct , the bound on instability index from Theorem |3.4| and the computed distribution of all 
purely imaginary gap eigenvalues. The figures in parentheses indicate the number of times a count is 
repeated, so that "8(3x), 7(6x)" in the last row of the table means that there are 8 gap eigenvalues 
for each j — 0, 1, 2, and 7 eigenvalues for each j = 3, 4, 5, 6, 7, 8. Accordingly, the total number of 
gap eigenvalues for m — 1,2,3,4,5 respectively is 22, 44, 78, 154, 268. The bound computed from 



Lemma 3.5 on the real part of any unstable eigenvalue satisfies A max < 1 in all cases. The value 
A = is always found to be a double zero of Ej for j '• = and 1, which corresponds to the four zero 
modes enumerated in Appendix B. 

For each m = 1, . . . , 5 the mechanism for transition to instability as uj decreases below uj cr is the 
same. A pair of purely imaginary gap eigenvalues, always with index j = 2, collides when uj = uj cr 
at A = A cr (as shown in Table 1) and move away from the imaginary axis in opposite directions, 
creating an unstable eigenvalue. Figure || shows how the graph of the real- valued function E%{it) 
(—uj < t < uj) changes as uj passes through the critical frequency for the case of spin m = 1. Figure ^ 
shows the critical wave profiles for m = 1, . . . , 5. 

Figure |l0| shows a log-log plot of uj cr versus m. (The transition frequency is calculated for 
noninteger values of m by solving the profile equation ( |l.4[ ) and computing the associated Evans 
functions.) From this data it appears that 

- 18 

uj* - w cr ~ — T (5.5) 

as m becomes large, suggesting that stable spinning waves may exist for any spin m. 

After this paper was submitted, the existence of a stability transition frequency uj ci for m = 1 and 
2 was also reported in a paper of Towers et al. ]3l| , who solve numerically the algebraic eigenvalue 
problem arising from discretization of a system of equations equivalent to ( |3.15 ). These authors also 



report finding (at the limit of their numerical accuracy) a weak instability for j = ±1 (see Figure 5 
in |3l]]) which affects the location of uj cr for m = 1 and 2 and causes them to find all waves with 
m = 3 unstable. If our present numerical results are correct, this weak instability is likely to be a 
numerical artifact, since for j = 1 the existence of a double eigenvalue at A = is consistent with 
the analysis of zero modes in Appendix H, and we find no other zeros of E\ (A) inside or near the 
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contour 7 except on the imaginary axis, where Ej(\) is real and changes of sign clearly locate zeros. 

We examined a few cases involving nonlincaritics other than cubic and cubic-quintic to see 
whether spectrally stable localized standing waves may exist. For the focusing-defocusing nonlin- 
earity g(w) = w 3 — w 7 , the critical frequency at which F has two zero-height hilltops is to* w 0.2722. 
Our computations indicate that the wave with spin m = 1 for uj = 0.24 is spectrally stable (j max = 9 
and |ReA| < 0.91 for unstable eigenvalues in this case). For the focusing-defocusing nonlinearity 
g(w) = w 5 — w 7 , the critical frequency « 0.0878, and there is a spectrally stable wave with spin 
m = 1 for u = 0.07 (j max = 11 and | Re A| < 1.27). 

For the saturable nonlinearity g(w) = w 3 /(l + w 2 ), which is always focusing, we found no 
spectrally stable spinning waves. Spinning solitary waves exist for < u < 1, growing in amplitude 
without bound as u) increases. With spin m = 1, we find for small lu that there are three unstable 
eigenvalues, as in the case of the focusing cubic above. For uj > 0.35, however, we find only one 
unstable eigenvalue; it has twist index j = 2. And as u> increases towards 1, the real part of the 
unstable eigenvalue becomes small. E.g., for u = 0.5 the unstable eigenvalue is A = 0.174 + 0.199i, 
and for ui = 0.9 it is 0.0372 + 0.0547i. As uj increases, then, we see that the nonlinearity saturates 
and the wave becomes less unstable. 



6 Discussion 

Our computations support the conjecture that for any focusing-defocusing nonlinearity, solitary 
waves with spin (encapsulated vortices) are spectrally stable for any value of the spin index m if 
they have no nodes (n = 0) and if the standing wave frequency w is sufficiently close to the kink 
frequency which is determined by the property that the graphs of y = g{x) and y = lu*x enclose 
two regions of equal area for < x < a*. (Equivalently, the potential F = F(s) has two zero-height 
maxima, at s = and s = a*.) 

To begin to understand why these waves might be stable, we observe that in the regime when 
uj ~ lu* and \m\ is large, the wave profiles develop a structure that may be accessible to analysis. 
This structure is suggested from Figures |] and 0, and a formal analysis is given in Appendix A. In 
this regime, the wave amplitude w(r) is approximately zero in a large central core region, then rises 
quickly in a transition layer where r ~ Ri n to a wide region where w(r) w a*, then decreases rapidly 
in a transition layer where r ~ i? ut to zero. 

Since the wave amplitude remains approximately constant near a* in a wide annular region 
where r is large and 9 is slowly varying, the solution is slowly varying in space. Now because the 
nonlinearity is defocusing for a mp litudes near a* (h! < 0) the spatially constant solution u = a^e 1 "** 



is spectrally stable, by remark 3.2. The zero solution is also spectrally stable, but spatially constant 
solutions with amplitude in the focusing range < a < a are not spectrally stable. 

In the zones of rapid transition where r is large and 6 slowly varying, the wave profiles w(r) 
appear to be well-approximated by a one-dimensional kink solution of the nonlinear Schrodinger 
equation ([hi]). Such a solution has the form je lul * t (l)*(x), where 0* is as described in section 2 and 
where 7 is a complex constant of unit modulus. This kink connects the zero solution to the solution 
7a >t e l "* t . It is significant to n ote that the kink solution itself is a spectrally stable solution of (fLl|). 
Because of the invariance of (1.1), we may take 7 = 1 in analyzing the stability of ^e w * *</)*. 



Proposition 6.1 Assume that g(u) — h(\u\ 2 )u where h is C 1 with h(0) — 0, and that properties 
(i) and (ii) of the introduction hold. Then admits a one- dimensional kink solution u(x,y,t) 



e 



lu} * t 4>*(x), where </)* is the unique solution of (\2. 4) and (2.t ), and this solution is spectrally stable 



Proof The existence and uniqueness of <fr„ are easy to establish. To analyze the stability, we perform 
the linearization of section 3 for uq(x, y,t) — e lul * t (/)*(x), obtaining dt'i' — A^>, where now 



A = J (A - a 2 + a(<f>*(x)) + /3(cj>*(x))R) 



(6.1) 
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Taking the Fourier transform in the transverse variable y, we obtain dt^ = A^S with A = —J(L+£ 2 I), 
where £ is the Fourier transform variable and 

T = (L+ V\ = (-dl-m*) 

l0 Lj V -dl-U{^)/^ 

Consider the spectrum of L. We know /*(0) < 0, /*(a*) < 0, and /*(a*) = 0, therefore 

lim /^(x)) < and lim f*{<t>*{x)) / 4>*(x) < 0. 

x — >±oo a: — *±oo 

Hence the essential spectra of L + and L_ are contained in [0, oo). Each operator L + and L_ 
annihilates a function of definite sign: 

L_ (f>*= and = 0, 

so a standard argument based on the maximum principle implies that L + and L_ have no negative 
eigenvalues. Thus L + £ 2 I has no negative eigenvalues, and it is simple to show that then J{L + £ 2 ) 
has no eigenvalues with nonzero real part. It follows that the spectrum of A is purely imaginary. 
Thus the kink e i "**</>„(a;) is spectrally stable, as claimed. D 

The focusing-defocusing properties of the nonlincarity thus permit the existence of a spectrally stable 
one-dimensional kink. 

A simplified description of the wave amplitude profiles, for spinning solitary waves when u) ~ lj* 
and \m\ is large, is that they have two pieces: 

(i) an outer kink region where w(r) 4>*(r — i? ut), 

(ii) an inner kink region where w(r) ~ <f>*{R m — r). 

Both the inner and outer kink radii R m and R on t become large, with i? ou t much larger that i?i n - 
Since both pieces approximate a spectrally stable one-dimensional kink solution, it becomes plausible 
that the whole wave can be stable. 




A Asymptotic analysis of wave structure 

As indicated in section 2, here we formally analyze the structure of nodeless localized waves in two 
limiting regimes, namely (i) when the spin m becomes arbitrarily large, and (ii) when the nonlinearity 
is of "hilltop" type and the standing wave frequency uj approaches the critical value lo* at which the 
positive maximum value of the potential F becomes zero, i.e., becomes equal to the local maximum 
at 0. 



A.l Nodeless wave location for large spin index m 

We observe numerically that the nodeless standing wave profiles w for large values of m appear to be 
copies of a single profile, shifted further to the right as m increases. To understand and characterize 
this property, we note that r is large where the wave amplitude is non-negligible. The term w'/r 
in (1.4) is thus apparently negligible, and we surmise that the quantity m 2 /r 2 is approximately 
constant where the wave amplitude is non-negligible. 

We therefore conjecture that w(r) ~ <pA (r — R) , where tpA is a putative positive solution to 

0"-A0 + /(0) = O (A.l) 



with maximum at the origin, that vanishes as r — > ±oo, and where A — m 2 /R 2 , with R to be 
determined. Since A > and /'(0) < 0, the solution cpA, if it exists, has exponential decay as 
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\r\ — > oo . Energy considerations show that such a solution <p A will exist provided that the even 
function Fa{s) = F(s) — \As 2 has a positive zero. 

To determine R or equivalently A, we derive an identity valid for profiles. Multiplying the 
equation ( p_.4| ) by w' and rearranging terms, we obtain 



1 r 
-W 

r 



1 r 

2 W ' 



—w 2 +F(w) 



= 0. 



Because w(r) ~ e?r m as r — > + , if m > 1 the quantity in square brackets vanishes as r 
it also vanishes as r — > oo. Integrating QA.2| ), we find 

™2 



f 



-w ' 



dr = Q 



for m > 1. 



Approximating r in the integrand of (A.3) by R, we find 



.4 



(f> A (r - R) 2 dr : 



4>' A {r - R) 2 dr. 



Extending the integration to the whole line, we arrive at the approximation 



.4 



4>a{x) 2 dx 



4>' A (x) 2 dx. 



(A.2) 
0+, and 

(A.3) 

(A.4) 
(A.5) 



We may compute the values of these integrals in terms of A by quadrature, as follows. Assume 
that Fa has a positive zero, and let B be the smallest positive z ero o f the function Fa- (That 
is, B is the solution to F(B) = ^AB 2 .) If <fi(r) is a solution to (A^l), then the energy E[4>] = 
^<fi'(r) 2 + FA{(j){r)) is independent of r. It follows from the facts that Fa(0) = and 4>A{r) — > 
as r — > oo that = and hence that the turning point of the solution (f>A at r — occurs at 

amplitude <\>a (0) = B. Furthermore, we may solve E = for <fi' to obtain 



|^| = y/-2F A (<f> A ). 

By changing variables from r to p = ^(r), wc find that 



/"Co /*.B 

' 0^(r) 2 dr= / ^-2F A (p)dp, 
o Jo 



dp. 



(A.6) 



(A.7) 



(A. 



y/-2F A (p) 

Since <j) A is an even function we may double these valu es to obtain the values of the integrals over 
the whole line. The resulting expressions inserted into ( A.5 ) yield an equation that determines A. 

As an example, we consider the cubic nonlinearity F(s) = —s + s 3 . We have F A (<fi) = jcj) 4 — 

+ A)<j>, hence B = ^2(1 + A) and 



and 



cf> A {r) 2 dr 



V2(T+A) 



^Ja+A)p 2 



hp 4 dp 



(l + Af 2 



V2(l+A) 



kP 4 



dp = 2Vl + A. 



(A.9) 



(A.10) 



Therefore, in this example, equation (A.5) yields A = -j, and we predict that the large-m profile will 
be centered at R — \/2 m and will have amplitude B = V3- This agree s wel l with the data shown in 
Fig. H We remark that in this example the values of the integrals in (A^) can also be determined 
from the explicit formula 4>a(x) = a/2(1 + ^4) sech(-\/l + Ax). 
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A. 2 Kink location in the flattop regime 

Here we st udy the structure of localized waves for "hilltop" potentials such as arise for the cubic- 
quintic in ( |l.2| ), when u> approaches the critical value uj* from below. For uj near w*, in the cubic- 
quintic case we observe numerically that the nodeless wave profile w(r) remains near w max ~ a* in 
a large interval, then decreases rapidly to zero when r is near the last point R where w(R) = \a*- 
Near this point the profile w(r) seems well- approximated by <f>*(r — R), where 4>* is the decreasing 



kink solution to ( |2.7| ) that satisfies (2 



To approximately determine R, we return to the identity (A. 2) and integrate from R — X to 
R + X to obtain 



R+X 



R-X 



1 



— W 

r 



— —w I dr = — 



1 ,n m 2 o . . 

2 W ^ W +F{W) 



R+X 



R-X 



(A.ll) 



We fix X large compared to 1 but independent of R. 



We first approximate the left side of (A.ll). Because w'{r) is nonnegligible only for r ~ R, and 



because w(r) « for r > R and w(r) « a* for r < R, we have 



R+X 



R-X 



-W ' 



w dr 



R+X 



w' 2 dr — m 2 (a») 2 



R-X 



1 

R-X r3 



■ dr. 



(A.12) 



We now approximate w'(r) in the integrand by cf>J (r — R) , and extend the integration to the whole 
line. Setting 

poo 

4>J{xfdx -- 



T 



we find 



R+x 



R-X 



-IV 



12 



in 



w 2 ] dr 



T 1 

R ~ 2 



y/-2F,(p) dp, 



m 2 (a«) 2 



(R-X) 2 



1 

R 2 



We next approximate the right side of (A.ll). Because w(r) w for r > R and w(r) 
r < R, we have 



1 



^w 2 + F(w) 



R+X 
R-X 



F ( a *) 



rrr 



2 (R-X) 2 



(a,) 2 



(A.13) 

(A.14) 

'• a* for 

(A.15) 



r(a*) 2 



2 (i? — A) 

where we have used t he fac ts tha t F*(a*) = and — F*(w) = ^(w* — w)u> 2 . 

Combining ( A.14 ) and ( A.15 ), we obtain the quadratic equation 



(w* - lo)R 2 - 



2T 

(a*) : 



-R — 777 = 



for i?, which depends only on uj and the nonlinearity /* . The solution for R is given by 

R ~ -Rout 



T+^JT 2 + (a*) 4 m 2 (uj* - uj) 



(a*) 2 (uj* — cj) 



(A.16) 



(A.17) 



In a similar way, we may approximate the location R of the inner kink of the nodeless wave 
profile with m > 2 for uj « lj* . We find 



-T+ y/T 2 + (fl,) 4 m 2 K -cJ) 
(a st ) 2 (cj st — uj) 



(A.18) 
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Note that 

2T 

Rout — Rin = 7 — 7^7 r, (A. 19) 

(a*)' ! {u!* — lu) 

independent of m. For fixed m > 2 we have 

i?i„ -> V ; and i? out -» oo (A.20) 

as lu approaches lu* from below. 

As an example, we consider the cubic-quintic nonlinearity f(s) = —lus + s 3 — s 5 . We have 
F(s) = -\lus 2 + \s A - is 6 , lu* = ^, a » = V3/2, and 

n/3/2 , o /o 

64 

Thus we obtain 



T= / ^2K{p)dp=\ J\jfi- l^ + ^dp^-^. (A.21) 



_ V3 +v /3 + 256 (3/16 - m 2 
^ out ~ 16(3/16-a;) ' [ ' 



-V3 + v /3 + 256 (3/16 - c) m 2 
^ = 16 (3/16 -c) ' (A - 23) 

which agrees well with the kink locations for the numerically computed wave profiles. E.g., for 
m = 3, ui = 0.18 we get i? Q ut ~ 52, i?; n ss 23 In the limit lu —t lu* note that R m — > 8m 2 / \/3. 



B Zero modes 

Solitary waves are not stable in the ordinary sense that the solution to the evolution equation with 
perturbed initial data differs little from the unperturbed solution for all time. This can be seen 
by considering, for example, two solitary waves with slightly different velocities whose initial data 
are approximately equal but which later drift apart. Instead, solitary-wave stability is defined in 
terms of the difference between perturbed and unperturbed solutions modulo symmetries of the 
evolution equation. For instance, the two solitary waves just discussed do remain close modulo 
spatial translation. 

Invariance of an evolution equation under symmetries gives rise to an eigenvalue-zero eigenspace 
for the generator of time evolution for the equation linearized about a solution. Linear instability 
due to generalized eigenfunctions in such a space need not correspond to a real instability, but rather 
to "drift" that can be cancelled by a symmetry transformation. On the other hand, an imaginary 
eigenvalue with a generalized eigenfunction (of index larger than 1) that does not correspond to a 
symmetry operation would suggest instability. 



In this section we make use of symmetries to enumerate the zero modes of (1.1), obtaining 
generalized eigenfunctions of the operator A corresponding to eigenvalue A = 0. We identify the 
same number of zero modes that we find numerically. 



Given a family u T of solutions to (1.1) with smooth dependence on parameter r, we can differ- 
entiate the identity 

— idtu T — Au T = g(u T ) (B-1) 
with respect to r to obtain a solution u to the corresponding equation linearized about Uq: 

- id t u -Au= a(\u \) u + (3(\u \)^-u (B.2) 

where 

u = d T u T \ T=Q . (B.3) 
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We apply this procedure to various one-parameter families of solutions, each having uq = 
e wt Vo(x) = e %u}t e lm0 w(r). With v given by u — e iult v, equation ( |B.2| ) becomes 

d t v = Bv (B.4) 

where B is the real-linear operator given by 

Bv = i ((A + a(w) -a 2 )v + (3{w)e 2im6 v) . (B.5) 



We note that (B.4) is equivalent to equation ( |3.4[ ) and that B is the single-component form of the 
operator A. (Because we are considering the real eigenvalue A = 0, in this section we prefer not to 
work with the two-component form and subsequent complexification.) 

• Invariance of NLS under phase changes: 

We define u T = e iT u Q . Because g(e lT u ) = e lT g(u Q ) and u is a solution of ([iTl|), u t is a 
solution for all real r. For this family, u — iuo and v — wq. Since v is time-independent, 



we see from (B.4) that v is an eigenfunction of B with eigenvalue zero. The corresponding 



eigenfunction = (Rev,Imv) of A lies in the j — subspace, according to (3J5), ( 3.13 ), and 



(P-14D 



• Invariance of NLS under rotation: 

We set u T = R t uq, where R T is the coordinate transformation that rotates uq by angle r. 
Thus u T = e luJt e lm< - e ~ T }w(r), whence u = —imuo and v — —imvo- This eigenfunction v of B 
is a real constant multiple of that induced by invariance under phase change. 

• Invariance of NLS under time translation: 

We set u T {x,t) = uq(x, t — r). Then u = —iluuq and v = —iluvq. This eigenfunction v of B is 
a real constant multiple of those induced by invariance under phase change and rotations. 

• Invariance of NLS under spatial translation: 

Writing (x, y) for the components of the spatial coordinate vector, we begin by setting u T (x, y, t) 



uo(x — T,y,t). Then u — 8 x uq and v — d x VQ. Since v is time-independent, we see from (B.4) 
that v is an eigenfunction of B with eigenvalue zero. By explicit computation, we see that 
according to (3.9), ( 3.13| ), and (3.14), the corresponding eigenfunction of A has components 



in the j ' = 1 and j = — 1 subspaces. 

In a similar way, setting u T (x,y,t) = uo(x,y — r,t) yields the eigenfunction v — d y vo of the 
operator B with eigenvalue zero, and the corresponding eigenfunction of A has components in 
the j = 1 and j = — 1 subspaces. It is straightforward to check that these two eigenfunctions 
resulting from translations along the x and y axes are linearly independent. 

Invariance of NLS under change of standing-wave frequency: 

As mentioned in Section ||, localized standing-wave solutions uq — e luJt e lme w(r) of (1.1) exist 



for a range of frequencies u. The radial profile w(r) depends on uj through the differential 
equation ( |l.4] ), in which f(w) — g{w) — low. Consequently, vq as well as uq is parametrized by 
to. Differentiating with respect to parameter w, we have 

u = itu + e lu) t d UJ v , v = itv + d^v . (B.6) 



Equation (B.4) thus gives Bv = ivo- Since we established earlier that B(ivo) — 0, we see that 



B (d w va) — ivq. Thus (again since B(ivo) = 0), 

B 2 (d^vo) = 0. (B.7) 
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Table 2: Zero modes 



Transformation 


\j\ 


Index 


Generalized eigenfunction 


Phase, rotation, time translation 





1 


iv 


Standing- wave frequency change 





2 




Space translation in x direction 


1 


1 


d x v 


Space translation in y direction 


1 


1 


dyv 


Boost in x direction 


1 


2 


ixvo 


Boost in y direction 


1 


2 


iyvo 



That is, dujVo is a generalized eigenfunction of B of index 2 corresponding to eigenvalue zero. 
Because 

d u v =e ime d u w(r), (B.8) 

the angular dependence of the generalized eigenfunction d^v® is identical to that generated by 
invariance under phase changes, and as a result the corresponding generalized eigenfunction 
^ of A lies in the j = subspace. 

Invariance of NLS under velocity boost: 

Let c € R 2 be a constant vector. If u(x, t) is a solution to then so is the function 

U c {x,t) = cxpji (\c-x- ^l c | 2 ^| u(x-ct, t). (B.9) 

Again writing (x, y) for the components of the spatial coordinate vector, we first let u T be the 
standing wave u$ boosted with speed r in the x coordinate direction (that is, c = (r, 0)): 

u T (x,y,t) = exp<i(^Tx-^T 2 tj \ u (x-rt, y, t). (B.10) 



Then u — ^xuq — td x iio and v = ^xvq — td x vo- Equation (B.4) thus gives Bv = —8 x vq. Since 
we established earlier that B (8 x Vq) = 0, we see that B(ixvo) = —2d x vo. Therefore (again 
since B (d x vo) = 0), 

B 2 (ixv ) = 0. (B.ll) 

That is, ixvq is a generalized eigenfunction of B of index 2 corresponding to eigenvalue zero. 
By explicit computation, we see that the corresponding eigenfunction ^ of A has components 
in the j = 1 and j = — 1 subspaces. 

In a similar way, letting u T be the standing wave uq boosted with speed r in the y coordinate 
direction (that is, c = (0,t)), we find that iyvo is a generalized eigenfunction of B of index 2 
corresponding to eigenvalue zero, and the corresponding eigenfunction of A has components in 
the j = 1 and j = — 1 subspaces. The two generalized eigenfunctions ixvo and iyvg resulting 
from velocity boosts in the x and y directions are linearly independent. 

We have thus found two j = generalized eigenfunctions, and four \j\ = 1 generalized eigenfunc- 
tions, for eigenvalue A = 0, as summarized in Table 2. 

To establish that these Jordan chains of generalized eigenfunctions do in fact terminate with 
index 2 in the general case, we make use of the following observation. 

Lemma B.l Suppose A = JL where L is self-adjoint and J* = J -1 = — J. Suppose ^ k is a 
generalized eigenfunction of A with index K corresponding to eigenvalue zero. If there exists $ in the 
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domain of A such that A$ = & K , then | J^ K+1 ^j ) = for j = 1, . . . , K, where &j = A K ~ j *S> K . 

Proof From the definition of ®j it follows that A^j = and A j ^fj = but ^ 0. Thus 

tyj is a generalized eigenfunction of A with index j . 

Define fy* = J^j- From the assumed form for A we then have 

A*tyj = —JAtyj = -J^j-i = (B.12) 

It follows that 

{A*yy* = (-lyjA*! = o. (B.i3) 

Thus, for 1 < j < K, we have 

(*, |J*K+l-i > = - \*K+l-i ) = - <** \A^ K ) 

= - (** lAi-^AQ) = - (** (B.14) 
= -((A*)J'*5|$)=0, 

as asserted. □ 

We now apply this observation to our index-2 generalized eigenfunctions. Consider first the 
generalized eigenfunction -02 = ixvo of B arising from velocity boost in the x direction. Then 
ipi = B?p2 = — 2d x vo and we have 

(*i |J# a ) = ^ <0i |t fa) = Re <-29^ |j z^ ) = - ||w || V 0. (B.15) 

Thus there is no generalized eigenfunction in this Jordan chain with index 3 (or higher). Similarly, 
analyzing the index-2 generalized eigenfunction ip2 = iy v o of B arising from velocity boost in the y 
direction, we find that there is no generalized eigenfunction in this chain with higher index. 

We next consider the index-2 generalized eigenfunction ip2 = d u v$ of B arising from change of 
standing- wave frequency which lies in the j = subspace. We set tpi = Bip 2 = ivo, and we find 

(*i |J*a> = Re(V»i \i*h) = Re(iv \i d u v ) = \d w \\v Q f . (B.16) 

Generically, we expect the quantity d w ||wo|| 2 to be nonzero, hence that i[>2 = d^VQ is the highest- 
index generalized eigenfunction in this Jordan chain. Our numerical studies confirm this expectation 
for nonlinearities other than pure cubic. 

The six-dimensional eigenspace for eigenvalue zero summarized in the table is in agreement with 
our numerical observations. For nonlinearities other than pure cubic, we find numerically that the 
eigenvalue zero has multiplicity two in each of the j = 0, j = 1, and j = — 1 subspaces. 

In contrast, in the special case of a purely cubic nonlinearity g(u) = 7|it| 2 u in two spatial 
dimensions, the norm ||uo|| is independent of u). To see this, we note that the function -Do defined 
by vq{x) = y/ui/j vq(t/u) x) satisfies the ^-independent equation A-Do — vq + |wo| 2- 5o = 0. It is 
straightforward to compute that, in N spatial dimensions, ||wo|| 2 = -cu^ 2 ~ N ^ 2 ||wo|| 2 . Thus, for the 
pure-cubic nonlinearity, in N = 2 spatial dimensions, the index-2 generalized eigenfunction does not 
necessarily terminate the Jordan chain in the j — subspace, and in fact generalized eigenfunctions 
with indices 3 and 4 are present in this special case. 



The structure of the generalized null spaces of A in the case of (1.1) with pure-power nonlinearities 
is established in Appendix B of the paper |33|. The index-3 generalized eigenfunction ip3 — ir 2 VQ 
of B arises from the invariance of cubic NLS in two spatial dimensions under the Talanov lens 
transformation. An index-4 generalized eigenfunction i/j^ has the form ip4 = £(r)e 4m61 , where £ is the 
bounded solution to the ordinary differential equation = r 2 WQ, where 

d 2 1 d m 2 2 
Ls = q^2 + ~ 72" ~~ u + 3 ^ w o- 

There are no generalized eigenfunctions with higher index. Our numerical observations in the pure 
cubic case show that eigenvalue zero has multiplicity four in the j = subspace, in agreement with 
this analysis. 



C Two proofs 

It remains to prove two technical results needed to characterize eigenvalues. 
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C.l Proof of Lemma [4.1 



To prove Lemma 4.1, we will transform the system (4.6) 

y' = B(r,j, X)y 

o _ o _ 

to an equivalent system of the form = (C +C)C where C is constant and C is integrable, then 
make use of facts concerning asymptotic behavior of solutions to such systems. To do so, we make 
the change of variables C±( x ) — h(r)y±(r) 1 where x ~ x(r) and h(r) are smooth monotonic functions 
such that 



x{r) 



lnr 

r 



for r < i, 



h(r) = 



for r > 2, 

With this change of variables, system ( |4.6| ) is equivalent to 

C = c{x)c, 



1 for r < §, 
Jr for r > 2. 



where 



C(x) = 




and 

Note that 



for r < | and r > 2. 



= M (r)y(r), M(r) = 



C(x) = r x {MBM- x + M'M~ 




h 
h'r x 







1 
r 



resp. 



l 

1s/r 



(C.l) 
(C.2) 

(C.3) 

(C.4) 
(C.5) 



We can write C(x) = C +C(x), where 



C = 




(J — my 



resp. 




for x < — ln2 and x > 2, and where the entries of C(x) are 0(l/x 2 ) as x 
x — * — oo, so that C(x) (but not xC(x)) is absolutely integrable on R. 

We discuss the asymptotic behavior of solutions as \x\ — > oo in the intervals (2, oo) and (— oo, — In 2) 

o 

separately. We note first that because fc+ and k- are nonzero, the matrix C is constant and diag- 

onalizable on (2, oo), with eigenvalues fc_ 
the fundamental set of solutions 



-k+, and — fc_. The system ^' = CC on (2,oo) has 



1 



V o 






(C.7) 
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It follows from Theorem IV. 1 of that for each member of this fundamental set, there exists a 
solution £ of ( |C.2| ) that is asymptotic to it as x — ► oo. Undoing the change of variables on the 
interval (2, oo), we see that system (fOj) has corresponding solutions respectively asymptotic to 



1 



o 
V o 



/ 



1 

Vft_ 





as r — > oo. Since the modified Bessel functions I v and K u have the asymptotic forms 



/„(*) 



1 



'2-kz 



1 



(C.8) 



(C.9) 



as \z\ — > oo in the right half complex plane, we see that the system (4.6) has solutions asymptotic 
to y i as r — > oo, for each i = 1, . . . , 4 as claimed. 

To analyze the asymptotic behavior of solutions as r — > + , we first suppose that both j + 777 

o 

and j — m are nonzero. Then the matrix C is constant and diagonalizablc on (— oo, — In 2), with 

o 

eigenvalues \ j + m\, \j — m\, — |j +m|, and —\j — m\. The system £ = C C on (— oo, — In 2) has the 
fundamental set of solutions 

/ 1 




,\j-m\x 




,-\j+m\x | -Ii+ m l 





,-\j-m\x 




(CIO) 



It follows as before that for each of these, there is a solution £ of ( p.2[ ) asymptotic to it as x — > — oo. 
Undoing the change of variables on the interval (— oo, — In 2), we see that the system (4.6) has 
corresponding solutions respectively asymptotic to 



r 



\j+m\ 



|j+m| 

r 

o 

V o 



Ai-m\ 




„-\j+m\ 



1 

\J+rn\_ 

r 






„-|j-m| 




(C.ll) 



as r — > + . Since for v > the modified Bessel functions /„ and K v have the asymptotic forms 



;(1 + 0(N)), ^W = |^(l + 0(|z|)) 



r^ + i) 



(C.12) 



as z — * 0, we see that the system (4.6) has a solution asymptotic to y i as r — » for each i = 1, . . . , 4, 
as claimed. 

o 

If either j + m = or j — m = 0, then on (— oo, — In 2) the matrix C is constant and has Jordan 
canonical form with a block of order 2 corresponding to the eigenvalue zero. Suppose for definitcness 
that j = to > 0. Then on (— oo, — In 2) we have 



C = 






1 





Q 


4m 2 




















1 















(C.13) 



and the system = C C has the fundamental set of solutions 






( 










-2m\ 














\ 


i r 





(C.14) 
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Since the entries of the matrix C(x) decay exponentially to zero as x — > — oo, xC(x) is absolutely 
integrable on (—00, — In 2), and we may apply Theorem IV. 4 of to establish that for each member 
of this set of solutions, there exists a solution £ of ( |C.2D that is asymptotic to it as x —* —00. Undoing 
the change of variables on the interval (—00, — In 2), we see that the system (4.6) has corresponding 
solutions respectively asymptotic to 



„2m 



as r — > 0. Since Iq(z) ~ 1 and Kq(z) ~ — 

O 

again has a solution asymptotic to y i as r 



( 1 \ 







( 


1 \ 




2m /r 1 








—2m/r 1 


C) 











h 


lnr 1 


V / 








/ 


\l/rJ 



(C.15) 



lnz as z — > 0, we see that the system (4.6) in this case 
— > for each i = 1, . . . , 4. This finishes the proof of the 



Lemma. IZ1 



C.2 Proof of Proposition [4.6 



It is evident that for any given A £ C„ the exterior systems associated with (4.6) and its adjoint 



(4.13) have the respective solutions y = yj ' Aj/f and z — z'f ^ a4°°' with the desired asymptotic 



behavior, since 



A y 2 0) 



A z 



(00) 



Vi a y 2 

o o 

z\ A 22 



as r 
as r 



(C.16) 



It remains to establish the uniqueness and analyticity of these solutions. Our plan is to transform 
the associated exterior systems to corresponding systems on the whole line by changing variables as 



in the proof of Lemma 4.1 Under this transformation, the quantities A and z^ 00 ' A z 



(00) 



map to solutions having asymptotic normalizations that enjoy a maximal rate of decay property. 
Uniqueness and analyticity of these sol ution s is a cons equence of the theory of |^6| . 

Let x(r), h(r), and M(r) be as in (C.l) and ( |C.3|) . Suppose y(r) is a solution of the exterior 
system 

y' = B A (r)y (C.17) 



associated with (4.6). Then y may be written as a linear combination of solutions of the form y Ay 
where y and y are solutions of (4.6). (Do this first for initial data, then solve the initial- value 
problem.) Given such solutions, then £(x) :— M(r)y(r) and £(x) := M(r)y(r) are solutions of 



(C.2), and C '■= C A C is a solution of the exterior system 



C A (x)C 



(C.18) 



associated with (C.2). Note that 



C A C = (My) A (My) = M(y A y), 

where the 6x6 matrix M has columns determined by M(ej A e^) = (Me/) A (Mejt). It follows that 
if we let C(x) — M(r)y(r), then £ is a solution of (C.18). Similarly, the converse holds: if ^ is a 
solution of ( |C.18| ) and y(r) = M{r)- l C{x), then y is a solution of QC . 1 7| ) . (Note that if JV = M _1 , 
then NM(y Ay) = y Ay for all 4-vectors y and y, hence N — M~\) 

We transform the adjoint exterior system in similar fashion. Suppose ff(x) = i;(r)M(r) _1 . Then 
z is a solution of the exterior system 



z' = -zB A (r) 



(C.19) 



associated with (4.13) if and only if f) is a solution of the exterior system 



V = -vC A (x) 



(C.20) 
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associated with the equation adjoint to (C.2). 

Let us now study solutions of (C.20) that correspond to solutions of ( p. 19 ) with the asymptotic 
behavior in ( C.16| ). Using the asymptotic relations established in Lemma 4.4, the definitions (4.15) 
and (4.4), and the asymptotics in (|C.9|), we determine that 



(z^ A z^M- 1 ~ (ziM" 1 ) A (z 2 M-^ ~ exp(-A (o °M™ (oo) 



(C.21) 



as x — > oo, where = fc+ + k— is an eigenvalue of Ca for x > 2, and ■u/ 00 -' is the particular 

associated left eigenvector 



«,(<»> = 



2^/k+k- 



(fef, 1,0,0) A (0,0, *_,!)• 



(C.22) 



Since fc + and fc_ have positive real parts for A e C CT , it follows that /i^ 00 ) is a simple eigenvalue, and 
^(oo) j g un iq Ue eigenvalue of largest real part. For A 6 C CT , it is clear that /}(°°) and fir 00 ) are 
analytic and the latter is nonvanishing. 

_ o , i 

Since Ca(s) = Ca (i)-Ca is integrable on K, we may now invoke the theory of |2q] . Proposition 
1.2 of |26| asserts (after space reversal) that ( C.20[ ) has a unique solution r) < - 00 - ) (a;, A) satisfying 



exp(fi ( - oc ^x)f] (oc) (x, A) -» ™ (oo) as a: 



CO, 



(C.23) 



and this solution is globally analyt ic in A for A 6 C ff . It follows that z'°°)(r, A) := ry^x , X)M (r) 
yields the unique solution of ( |C. 19 ) having the asymptotic behavior asserted in Proposition |4.6| , and 
that i^ 00 -* is analytic in A f or A S C CT . 

The treatment of ( C.17 ) is similar, at least in the case when |j ; + m| and \ j — m\ are both nonzero, 
which we consider first. As x — > — oo we have 



M(yf ] A yf ) ~ (M&) A My° 2 ~ exp(A ( -° o) x)i> 



(C.24) 



where /t*- -00 -* = |j + m| + |j — m| is an eigenvalue of C A for x < — In 2, and tr °°' is the particular 
associated right eigenvector 



r(|j + m|)r(|i-m|) 




(C.25) 



Clearly /t( 00 ' is independent of A, and fr is analytic and nonvanishing in the cut plane C CT . 
Provided \j + m\ and \j — m\ are nonzero, /t^ 00 ) is simple and is the unique eigenvalue of largest 

o , , 

real part for C A for x < — In 2. Proposition 1.2 of (26) implies that ( C.18 ) has a unique solution 
£ (x, A) satisfying 



- (-oo) 

tr ' as i 



-oo, 



(C.26) 



and this sol ution is globally analytic for A e C a . It follows that y^ := M 1 ^ is the unique 
solution of ( C.17 ) having the asymptotic behavior asserted in Proposition L6, and that y (0) is 
analytic in A for A G C a . 

It remains to consider the case when either \j + m\ or \j — m\ is zero. For definiteness, suppose 
j = m > 0. We still have (C.24), and jj 1 - -00 ^ is the unique eigenvalue of largest real part, but now it 
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is a double eigenvalue of geometric multiplicity one, with eigenspace spanned by tr °°\ The proof 

of Proposition 1.2 of pfj must be modified to obtain the uniqueness and analyticity of £ 
This can be done as follows. First, we change variables in ( p. 18 ), writing 



(-oo) 



v{x) = exp(-A ( ~°°MC(20 (C.27) 

o _ 

Let B = C A — and i?(x, A) = C A (x). Then the first equation in ( p. 18 ) is equivalent to 

v'(x) = (B + R{x,X))v. 
Note that there is a constant C\ > independent of A such that 

|| exp(Bas)|| < Ci(l + |a:|) 

for all x > 0. Also, the entries of i?(x, A) decay exponentially as x — > — oo, uniformly for A in any 
compact subset of C a . For any compact subset Qi C C CT , there exists xo sufficiently negative so that 



C\ sup 



(l + |s|)||i?(s,A)||ds< 1. 



We define an operator T = J-(X) on the space C((— oo,xo]) of bounded continuous functions on 

(-oo, x ] by 

jR>(x) = [ X eS x - s)B R(s,\)v(s)ds. 

Then it follows that for A G Ox, 



sup |^u(x)| < 8 sup \v(x) 



x<x 



x<x 



From here, the proof goes the same as the proof of Proposition 1.2 of 1 26| . The equation v = 
v + J-v yields a bicontinuous correspondence between bounded solutions of (C.27) on (— oo,xq], and 



those of v = Bv, all of which have the form v(x) — cv ( ca> for some c € C. The solution we seek 
has the representation 



(x) =cxp(jj(-° D '>x){I - F)- 1 *) 



-!,-,(-«) 



for x < xo, and is extended by solving ( p. 18 ). It is independent of xo and is obtained by iterations 

~( — oo) 

that are analytic in A and converge uniformly for x < xo, A € It follows that £ (x) is globally 
analytic for A € C CT . This finishes the proof of the Proposition. D 
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